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ABSTRACT 

Type III solar radio storms, observed at frequencies below ~ 16 MHz by space 
borne radio experiments, correspond to the quasi-continuous, bursty emission of 
electron beams onto open field lines above active regions. The mechanisms by 
which a storm can persist in some cases for more than a solar rotation whilst ex- 
hibiting considerable radio activity are poorly understood. To address this issue, 
the statistical properties of a type HI storm observed by the STEREO/WAVES 
radio experiment are presented, examining both the brightness distribution and 
(for the first time) the waiting-time distribution. Single power law behavior is 
observed in the number distribution as a function of brightness; the power law 
index is ~ 2.1 and is largely independent of frequency. The waiting-time dis- 
tribution is found to be consistent with a piecewise-constant Poisson process. 
This indicates that during the storm individual type HI bursts occur indepen- 
dently and suggests that the storm dynamics are consistent with avalanche type 
behavior in the underlying active region. 

Subject headings: Sun: activity — Sun: corona — Sun: radio radiation — meth- 
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Introduction 



Type III solar radio storms consi st of many individual typ e I II rad io bursts produced 
quasi-continuously over several days ( Fainberg &: Stone 1970a| jbl. Il971 ). Such storms can 
persist for long intervals - they have been observed for up to half a solar rotation, and 
in some cases over a whole rotation ( Fainberg fc Stone 1970al ). It h as been shown that 
they are associated with both active regions and metric type I storms (IKayser et al. 1119871 : 
Gopalswamv I l2004r) ■ The type III storm emiss ion is produced by electron beams on open 
field lines ( Bougeret. Fainberg fc Stone 1984a| jbl). and so type III storm activity may shed 
light on the dynamics of the open/closed field line boundaries associated with the underlying 
active region. The dynamics of type III storms are still not well understood, particularly the 
mechanisms by which they are able to persist for such long times whilst at the same time 
generating considerable numbers of radio bursts. Fur thermore, type III storms have also 
been observed as precursors to coronal mass ejections (IReiner et al. II2OOII ). and so a better 
understanding of their properties may help to elucidate the dynamics of active regions prior 
to eruptive events. 

One important statistical prope rty is the distribution of pow er S in the individual bursts. 
In their initial type III storm study, iFainberg fc Stone I (11970a[) noted that there were more 
faint bursts than bright ones. More recently, iMorioka et al. I (120071 ) studied the statistics 
of type III bursts using Geotail (24Hz - 800 kHz) and Akebono (20kHz - 5MHz) radio 
data. Using observations covering 4 Bartels rotations and containing both ordinary type 
HI emission and intervals of type HI storm activity, they found that the overall occurrence 
of type HI bursts followed a double power law; for example, at a frequency of 860 kHz 
and below an intensity of ~ 10^^^ V^m^Hz"^ (~ 2.7 x 10~^* Wm~^Hz"^), the number 
distribution (simply the number of bursts at a given brightne ss) followed a S~^'^ po wer law 



whereas above this limit, a 5""°'^^ power law was observed. iMorioka et al. I (120071 ) named 
the lower flux distribution 'micro-type HI bursts', and concluded that although they were 
associated with the same active regions as the 'ordinary' type HI bursts, they were produced 
by independent, coexisting processes. This is consistent with the association of the Type HI 
storms with meter- wave Type I 'noise storms', which do not have the close association with 
flares that characterizes ordinary type HI bursts (e.g. iGopalswamy II2004J ). In comparison. 



Fitzenreiter. Fainberg fc Bundy I (119761 ) found in a study of ordinary type HI bursts observed 
over a period of three months that the number/unit flux density fell off with increasing flux 
following a 5*"^ power law with (3 varying from 1. 3 3 at 1 10 kHz to 1.69 at 4.9 MHz. This in 
fact appears compatible with the IMorioka et al. I (120071 ) study, since a S"^'^^ power law in 
the number /unit flux density corresponds to a S""*^'^^ power law in the number distribution. 
For comparison, at higher frequencies (164 MHz, 237 MHz and 327 MHz), the number 
distribution of type I bursts in radio noise storms has been found to follow a similarly steep 



- 3 - 



S ^ power law (IMercier fc Trottet 1119971 ) 



A second statistical property is the Waiting-Time Distribution (WTD), that is the 
distribution of times between indivi dual type III events within the storm. WTD ana lysis 



Coronal Mass Ejections (jWheatland 



has been applied to soft X-ray flares (Wheatland II2OOOI : IWheatland fc Litvinenko II2OO2I ) and 



20031 ) and in both cases it has been shown that the data 



are consistent with a time-dependent Poission process, in which events occur independently. 
The WTD prov ides models with extra constraints; for example, flare avalanche models (e.g. 
Lu et al. II1993I ) predict that flares occur independently, and so the WTD is an important 
observable. 

To our knowledge the WTD has not been used to examine type III storms, particularly 
at the low radio frequencies typically measured by spacecraft (below ~ 16 MHz). In this 
Letter we analyze both the brightne ss distribution and the WTD of a type III storm observed 
by the S /WAVES ra dio experiment (IBougeret et al. II2OO8I ) on board the STEREO spacecraft 
(IKaiser et al. II2008I ). in order to examine in more detail the underlying processes responsible 
for type III storm emission. 



2. Observations 



Each STEREO spacecraft is equipped with an identical S/WAVES instrument designed 
to measure both in-situ plasma wav es and freely propaga ting radio emission in the fre- 
quency range 2.5 kHz - 16.025 M Hz jBoueeret et al. l EoOsf). S/WAVES uses three orthog- 



onal monopole electric antennas (iBale et al. II2008I ) connected to two instrument channels 



in a typical operating mode, one channel is connected to a pseudo dipole and the other to 
a monopole. Within each channel, there is a high frequency receiver (HFR) and a low fre- 
quency receiver (LFR). The observations presented here were obtained shortly after launch 
using the HFR which covers the frequency range 125 kHz - 16.025 MHz, before the lunar 
flybys that sent the STEREO spacecraft into their final heliocentric orbits and before many 
of the other instruments became operational. Consequently the radio data from the two 
spacecraft are essentially identical, and so here we present the analysis from STEREO A. 

Figure [1] shows measurements from the S/WAVES experiment on STEREO A measured 
between 0.2 MHz and 10 MHz from 8 November 2006 0000 UT - 17 November 2006 0000 UT; 

9 days of data are shown. The data are from channel 1 (corresponding to the - Ey dipole), 
and have be en calibrated in u r iits o f incident flux S [Wm~^Hz~^] following the procedure 
described by Eastwood et al. I (120091 ). The calibration procedure is most accurate below ~ 

10 MHz; at higher frequencies, the antenna resonance is less well modeled. Lines of Radio 
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Frequency Interference (RFI) at fixed frequencies have been removed by averaging across 
the frequency band. Terrestrial RFI is still observed periodically between 9.5 MHz and 10 
MHz; in the following analysis, frequencies have been chosen to avoid this interference. 

Type HI emission began to be observed intermittently as early as 9 November 2006, 
and the storm itself was established by 11 November 2006. Active region 10923 was the 
only a ctive region ori the disk during this time and so we associate it with the type HI 
storm (IKayser et al. 1119871 ) although it appears that the storm developed after the active 
region appeared on the disk. Over the next 6 days, the number, frequency and intensity of 
type HI emissions increased. On 17 November, the STEREO spacecraft executed an Earth 
perigee pass. Whilst the storm continued, S /WAVES observed overwhelming terrestrial RFI 
at all frequencies. Active region 10923 was close to the limb by this time, and a drop in 
the number and intensity of type HI bursts on 18 November was followed by the leading 
spot of AR 10923 rotating past the limb on 20 November. This confirrned the association of 
the active region with the type HI bursts (IBougeret. Fainberg fc Stone Ill984bl ). but suggests 
that the activity is more associated with the leading edge of the region. In this Letter the 
interval 11 November 2006 0000 UT - 17 November 2006 0000 UT has been analyzed. 



3. Analysis and Results 

To examine the properties of this storm interval in more detail, time series of data at 
specific frequencies were extracted and analyzed. Identification of individual type HI bursts 
was made in the following way. We first established a threshold brightness 5*0 and defined a 
burst as an interval during which time the brightness S satisfies S > Sq. The burst was then 
assigned a brightness S'burst = S'max ~ S^, the maximum brightness during the burst interval 



less the galaxy brightness Sg (ICane Ill979l ). The time of the burst is similarly defined as the 
time the maximum brightness was observed. Inspection of the data during the most active 
intervals shows that in general, bursts are typically observed to last only one timestep (~ 
25 s), and bursts lasting multiple time steps are rare. Accordingly mis-counting from event 
overlap is not significant. This results in a set of bursts and from their observed times, a 
set of A^ — 1 waiting times. 

We first consider the brightness distribution f{S), which we define as the number of 
bursts per unit brightness. If this distribution follows a power law 



f{S) = AS-' (1) 
then the maximum likelihood method (iBai Ill993l ) can be used to show that the most likely 
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value of the power law index 5m is given by 



N 



+ 1, 



(2) 



where Sj is the bri ghtness of the i^^ burst. The error associated with 5^ is given by 



flWheatland 112004 ): 



{8m - l)iV-V2. (3) 



Rather than examining the brightness distribution directly, it is preferable to plot the 
cumulative brightness distribution, which avoids binning of the data, and more easily shows 
the effect of the choice of 5*0. The cumulative distribution also follows a power law as shown 
by equation (jl]): 

j-OO A 

Ns:.So= AS-'dS = —-S,'^\ (4) 

This analysis was applied to the 3.025 MHz time series. Figure [2] shows the distribution 
of cumulative burst brightness (points). A threshold 5*0 = 5 x 10~^^ Wm~^Hz~^ is used, 
well above the average galactic brightness of ~ 1.1 x 10^^^ Wm~^Hz~^. This results in a set 
of 824 bursts. In this case, the power law index associated with the brightness distribution 
is (5 = 2.16 ± 0.04. The power law derived from the data using the maximum likelihood 
method is shown as a solid line in Figure El dashed lines show the power law plus/minus 
the error calculated using equation ([3]). The data in general follow the power law and the 
significance of any departure in the final few points is marginal; there is no substantial 
evidence for a break to a different power law index although a detailed quantitative model 
comparison has not been done. If a higher value of So is chosen, this has the effect of 
removing points on the left hand side of the cumulative distribution, but does not change 
the index of the power law. If a smaller value of 5*0 is chosen, the data eventually rolls over 
to the background level. This procedure was repeated for other frequencies and the results 
are shown in Table [H It can be seen that there is very little variation in the power law index 
as the freq uency changes. A larger power law ind ex is seen compared to the study of type 



III bursts (IFitzenreiter. Fainberg fc Bundy Ill976l ). 



The top panel of Figure [3] shows the time series at 3.025 MHz, where the horizontal 
line indicates the threshold used in the analysis. The middle panel of Figure [3] shows the 
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cumulative number of bursts as a function of time. It can be seen that the number of bursts 
does not increase at a constant rate, rather there are intervals where the number of bursts 
increases slowly, followed by shorter intervals where the burst rate is enhanced, and the 



( Scarele 


1998; 


Wheatland 


2000) 



consistent with a constant-rate Poisson process. Since a Poisson process has a waiting-time 
distribution given by 



P{At) = Xe 



-XAt 



(5) 



where At is the interval between events, and A is the rate, a piecewise constant Poisson 
process made up of M intervals with rates Aj and durations ti has a WTD given by: 



M 



P(At) 



1 \^ -A,At 

A Z^T^^^ 



(6) 



where A is the average burst ra te and T is the duration of the observing interval (e.g. 
Wheatland fc Litvinenko 1120021 ). The lower panel of Figure [3] shows the Bayesian blocks 
rate decomposition. The Bayesian procedure decided that there were M = 16 intervals 
with approximately constant rate during the observing period. 



The top panel of Figure H] shows the observed WTD as points. The data are again 
presented cumulatively to avoid binning of the data. The model cumulative WTD is given 
by: 



P{At > Ato 



P{At)dAt 



Ato 



:A,e" 



-XiAto 



(7) 



where the rates Ai and intervals ti come from the Bayesian blocks procedure. This model is 
also shown in the top panel of Figure HI The bottom panel of Figure H] shows the ratio of 
the data and the model as a function of waiting time. There is good agreement between 
the data and the model, which strongly suggests that the observed WTD is consistent with 
a piecewise constant Poisson process whose rate varies as shown in Figure [3l Similar results 
were obtained for the different frequencies shown in Table [H 
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Discussion 



As mentioned in the introduction, type III emission is thought to arise from non-thermal 
electron beams on open field lines. Emission may occur at the fundamental or its harmonic; 
type III bursts wh ose trajectory intersects the spacecraft show that the emission is at the 
fundamental mode Jheblanc. Dulk fc Bougeret lllQQsl: [Pulklboooh. but the possibility of har- 



monic emission perhaps should not be discounted (IKellogg 



IQSOuBougeret. Fainberg fc Stone 



1984bl ). Since the local electron plasma frequency depends on density, that is /p[Hz] 



8.98 X 10^ X (nc[cm ^])^/^, a density mo del is required to determine the h eight of the radio 



emission. Using the model developed by lLeblanc. Dulk fc Bougeret I (Il998l ). in this example 
2 MHz corresponds to ~ 3 Rs and 9 MHz corresponds to ~ 1.8 Rs- The model is scaled 
to the solar wind density at 1 AU; between 17-20 November 2006 the average solar wind 
density measured by the ACE spacecraft at 1 AU is approximately 4 cm~^. If the type HI 
emission is harmonic, then then the local density is reduced, and the emission comes from 
higher in the corona - between approximately 2.5 and 9 Rs for 9 MHz and 2 MHz. 

It is generally accepted that the radio emission arises from the non-linear conversion 
of Langmuir waves excited at the local plasma frequency by non-thermal electron beams. 
As the electrons move along the magnetic field, the faster electrons outrun the slowe r elec- 
trons leading to a bump-on-tail distribution and instability (IRobinson fc Cairns Il2000l ). The 
growth of Langmuir waves causes the distribution to relax to marginal stability; in stochastic 
growth theory, fiuctuations in the ambient density are thought to cause fiuctuations about 
marginal stability which allows wave growth without rapid destruction of the electron beam 
( IRobinson 1119921 : iRobinson. Cairns fc Gurnett Ill993l ). This suggests that the electron beams 
are being emitted from closer to the Sun. The more general association of type HI storms 
(open field) and type I storms (closed field) suggests that the storm activity is related to 
reconfiguration of the boundary between open and closed magn etic fields in the underlying 
active region. Analysis of Potential Field Source Surface maps (ISchrijver fc DeRosa II2003I ) 
shows that a small region of open field lines was persistently present adjacent and to the 
west of the active region during the storm. The interface between the open and closed field 
regions rotated behind the limb before the active region itself, which may help to explain 
the observed decline in storm activity prior to the active region rotating past the limb. 

In this type HI storm, a S^"^'^ power law was found in the number distribution, cor- 
responding to a S~^'^ power law in the number/unit fiux density. The power law index 
was found to be largely independent of frequency. Although the exact relationship between 
the primary energy release, the energy and properties of the electron beam, and the energy 
converted to radio emission is not well understood, it is likely that this r efiect s power-law 
behavior in the underlying dynamics. In a previous study, iMorioka et al. I (120071 ) observed a 
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^-3.6 power-law in the number distribution of the storm ('micro') burst population, but this 
was averaged over several storm intervals and used data binning. In this analysis, maximum 
likelihood was used, although this technique should be used with care if there are departures 
from the power-law behavior due to the method of event selection, background subtraction, 
etc. It is necessary to study other events to understand whether this difference represents 
variability from storm to storm, or is due to some other effect. The rate at which bursts are 
observed is found to change through the storm. Using a Bayesian block decomposition, we 
find the observed WTD to be consistent with a piecewise-constant Poisson process, which 
indicates that individual type III bursts occur independently of one another. It is interesting 
to note that the Bayesian block structure has no correlation with the occurrence of flares 
in the active region. A GOES C3.3 event occurred on 12 November 2009 at 2155 UT, for 
example, in the middle of a long Bayesian block. 

These two properties are predicted by avalanche models of driven systems. In the so- 
l ar context, such models have been developed to understand the properties of solar flares 



(ILu et al. 1 119931 ). This suggests that in the active region under consideration, magnetic en- 
ergy and flux, slowly injected from below and being driven e.g. by twisting due to foot point 
motion, is being released in a series of independent events which result in the emission of 
electron beams onto open field lines. In particular, the long lifetime of type III storms indi- 
cates that the active region can remain in such a dynamically balanced state for a prolonged 
period of time even as the burst rate slowly changes. This stability was not disturbed by the 
~ 21 GOES flares in this time interval from AR 10923, including 5 C-class events before 14 
November 2006. 
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Fig. 1.— STEREO- A S/WAVES observations of radio emission from 8 Nov 2006 OOOOUT to 
17 Nov 2006 OOOOUT (9 days). A type III storm associated with Active Region 10923 was 
observed, beginning on 11 Nov. Terrestrial RFI is observed periodically between 9.5 MHz 
and 10 MHz but is avoided in the analysis. The white line at 2 MHz is an artifact of the 
plotting software. 




Fig. 2. — Cumulative distribution of burst brightness at 3.025 MHz, together with the 
maximum likelihood fit. 
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2006-1 1-1 1/00:00:00 - 2006-1 1-1 7/00:00:00 
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November 2006, 3.025 MHz 

Fig. 3. — (Top) Time series at 3.025 MHz (middle panel), cumulative number of bursts 
(bottom panel), and results of the Bayesian blocks decomposition. 
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STEREO A Type III storm 
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Fig. 4. — The top panel shows the cumulative waiting time distribution of bursts at 3.025 
MHz, together with model based on a piecewise constant Poisson process, with rates deter- 
mined by the Bayesian blocks procedure. The bottom panel shows the ratio of the data and 
the model. 
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Table 1. Power law index as a function of frequency 



Frequency (MHz) 


Power Law Index 


3.025 


2.16 ± 0.04 


4.025 


2.13 ± 0.04 


5.025 


2.10 ± 0.05 


6.025 


2.15 ± 0.06 


7.025 


2.24 ± 0.06 


8.025 


2.22 ± 0.07 


9.025 


2.36 ± 0.07 


10.025 


2.20 ± 0.07 



